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ABSTRACT 

Stress- wave  profiles  in  vitreous  Ge(>2  induced  by  planer  and  spherical  projectile  impact 
were  measured  using  piezoresistance  gauges  in  the  4  to  1 8  GPa  shock  pressure  range.  The 
planer  experiments  demonstrate  that  the  response  of  vitreous  Ge02  to  planar  shock  loading  can 
be  divided  into  three  regimes:  (1)  An  elastic  shock  regime  with  ramp  Hugoniot  elastic  limit 
(HEL)  precursor  with  a  peak  stress  and  particle  velocity  of  4  GPa  and  0.33  km/s.  Its  propagation 
velocity  decreases  from  an  initial  longitudinal  elastic  wave  speed  of  3.5  km/s  at  essentially  zero- 
amplitude,  to  2.8  km/s  at  4  GPa;  (2)  A  transition  wave  regime.  Here  a  ramp  wave  is 
superimposed  on  the  precursor  with  an  amplitude  of  0  to  2  GPa  for  peak  loading  pressures  of  8  to 
14  GPa.  Above  4  GPa  the  ramp  wave  velocity  decreases  to  a  value  below  2.5  km/sec  (the  speed 
of  the  bulk  wave,  at  the  HEL);  (3)  A  shock  wave  achieving  the  final  shock  state  forms  when 
peak  pressure  is  >6  GPa.  This  wave  can  be  described  via  the  linear  shock-particle  velocity 
relation  D  =  0.917  +1.71  u  (km/s)  over  the  6-40  GPa  range  for  an  initial  density  of  3.655 

g/cm3.  Combining  the  Jackson  and  Ahrens  [Jackson  and  Ahrens,  1979]  data  and  the  present 
demonstrates  the  phase  change  from  4-  to  6-fold  coordination  of  Ge+4  with  O'2  in  Ge02  occurs  in 
the  pressure  range  of  4  to  15  ±  1  GPa  under  planar  shock  loading.  Comparison  of  the  shock  data 
for  fused  Si02  with  those  for  vitreous  Ge02  demonstrates  that  range  of  transformation  to  a 
denser,  presumably  6-fold  coordinated  structure  in  both  media  are  similar.  The  Hugoniots  of 
GeC>2  and  SiC>2,  both  initially  vitreous,  are  found  to  be  virtually  coincident  if  pressure  in  SiC>2  is 
calculated  by  multiplying  the  Ge02  pressure  by  the  ratio  of  the  initial  densities  of  vitreous  GeC>2 
to  fused  SiC>2.  The  volume  axes  are  translated  with  the  specific  volumes  upon  onset  and 
completion  of  the  phase  change  are  aligned.  This  result,  as  well  as  the  usual  crystal-chemical 
structural  arguments,  provides  the  basis  for  considering  vitreous  GeC>2  as  an  analogous  material 
to  fused  SiC>2.  Our  spherical  projectile  impact  results  demonstrate:  (1)  The  supported  elastic 
shock  in  fused  Si02  decays  less  rapidly  than  a  linear  elastic  wave  when  elastic  wave  stress 
amplitude  is  higher  than  4  GPa.  A  supported  elastic  precursor  in  vitreous  Ge02  decays  faster 
with  radius  than  a  linear  elastic  wave;  (2)  In  Ge02  (vitreous)  unsupported  shock  waves  decay 
with  peak  pressure  in  phase  transition  range  (4-15  GPa)  with  propagation  radius  (r)  as  f3'35, 

close  to  the  prediction  of  Chen  [1999].  Based  on  a  simple  analysis  of  spherical  wave 
propagation,  we  find  that  the  different  decay  rates  of  a  spherical  elastic  wave  in  fused  Si02  and 
vitreous  Ge02  are  predictable  on  the  basis  of  the  magnitude  and  sign  of  9K/9P  hug,  where  K  is 
bulk  modulus  and  P,  shock  pressure. 

INTRODUCTION 

In  the  1960’s  and  early  1970’s  when  the  first  shock  wave  measurements  were  conducted 
on  silicates  [Wackerle,  1962],  [Fowles,  1967],  [Ahrens  and  V.  G.  Gregson,  1964],  [McQueen  et 

al.,  1967],  [Ahrens  et  al.,  1969]  it  became  clear  that  this  class  of  materials  all  demonstrated 
behavior  which  indicated  that  transformation  to  dense  structures  in  which  Si+4  was  coordinated 
to  six  0‘2  ions  as  in,  the  then  recently  discovered,  (at  Meteor  Crater,  Arizona)  high-pressure 
phase  of  Si02,  stishovite  (rutile-structured  phase)  which  has  a  density  of  4.30g/cm3,  some  62% 

and  95%  greater  than  crystal  or  fused  quartz  starting  materials  [Chao,  1962].  Motivated  by 
scientific  curiosity  and  a  need  to  contain  explosive  products  and  understand  the  interaction  with 
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rock  of  nuclear  explosions  in  the  Earth,  a  detailed  series  of  measurements  of  release  isotropes  of 
SiC>2  and  other  silicate  minerals  shocked  into  the  mixed  phase  and  high  pressure  phase  regime 

were  obtained  [Grady  et  al.,  1974],  [Podurets  et  al.,  1988],  [Chhabildas  and  Miller,  1985],  All 
these  data  point  to  the  general  conclusion  that  transformation  to  the  high-pressure  phase  occurred 
upon  overdriving  the  low-pressure  phase  to  a  significantly  higher  shock  pressure  than  required 
for  thermodynamic  equilibrium  by,  in  several  cases,  a  factor  of  two.  Upon  unloading,  the  mixed 
phase  and/or  high-pressure  phase  regimes  demonstrated  hysteretic  behavior.  This  indicates  that 
unloading  was  occurring  in  pressure-density  plane  along  the  isentrope  of  the  high-pressure  phase 
and,  as  a  result  of  this  hysteresis,  the  silicates  appeared  to  be  extremely  absorptive  of  shock 
energy,  e.g.  Swegle  [1990].  Only  recently  have  we  understood  this  behavior  in  terms  of  the 

inclusion  theory  of  Truskinovsky  [1984],  which  although  developed  to  describe  quasi-static 
transformation  is  expected  to  be  applicable  to  many  shock-induced  transformations,  including 
the  quartz-stishovite  transformation  [Stolper  and  Ahrens,  1987]  because  it  is  martinsitic  in  nature 
as  are  other  major  and  important  phase  changes  such  as  the  a  -»  e  phase  (iron),  B 1  B2  phase 
(alkali  halides),  and  graphite  — >  diamond. 

Although  it  was  long  recognized  that  the  phase  transformation  in  silicates  and  SiC>2 
occurred  at  stress  levels  of  -10-20  GPa.  This  stress  level  and  hence  these  transformations  were 
out  of  range  of  interest  of  light  armor  materials  that  are  effective  against  small  arms  projectiles. 

However,  in  the  case  of  Ge02  and  other  germenates  this  does  not  seem  to  be  the  case. 

For  GeC>2  (which  has  a  phase  diagram  similar  to  Si02)  the  glass  (3.64  g/cm3)  to  rutile-like  phase 
(6.25  g/cm3)  transformation  demonstrates  a  density  increase  of  72%  that  occurs  over  the  pressure 
range  of  4-15  GPa.  This  shock  stress  range  is  inferred  on  the  basis  of  the  present  work  and  the 
early  studies  of  Jackson  and  Ahrens  [Jackson  and  Ahrens,  1979].  This  is  the  pressure  range  of 
interest  upon  small  arms  fire  striking  GeC>2  glass  and  hence  this,  in  combination  with  other 
materials,  could  possibly  be  useful  as  a  potential  ingredient  in  a  transparent  armor  system. 

In  this  report  we  briefly  review  the  previous  work  on  a  program  initiated  in  1996  in 
which  the  release  isotropes  of  Ge02  were  first  measured  and  these  were  used  to  model  excess 
attenuation  in  the  phase  transformation  region.  Our  recent  work  in  which  shock  structure  and 
attenuation  in  Ge02  in  one-dimensional  longitudinal  compression  and  a  spherical  diverging 
shock  wave  are  summarized. 


PREVIOUS  RESEARCH  ON  Ge02 

Motivated  by  the  possibility  that  further  coordination  increases  in  zM+402 ,  e.g. 

viM+402  v,im+4o2  -4  viiiM+402 

where  %  is  a  6-,  7-,  or  8-fold  coordination  of  a  metal  ion  (M)  in  column  4A  of  the  periodic  table 
(M  =  Si,  Ge,  Sn,  Pb)  or  in  column  4B  (M  =  Ti,  Zn,  Hf,  Rf),  Jackson  and  Ahrens  [1979]  studied 
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shock  compression  of  GeC>2  both  in  vitreous  (4-fold  coordinated)  and  in  rutile  (6-fold 
coordinated)  structures  to  determine  whether  further  transformation  occurred. 

The  column  4  elements  mentioned  above  actually  achieve  a  number  of  structures  and 
several  that  have  higher  coordination  than  6.  These  have  been  studied  and  classified  by  J.M. 

Leger  [1996].  The  Hugoniot  data  for  GeC>2  is  plotted  in  Figure  1  together  with  the  static 
compression  data  for  GeC>2,  starting  from  glass  and  quartz  structure.  Also,  GeC>2  in  the  quartz 
phase  transforms  to  a  vitreous  phase  upon  recovery  from  dynamic  compression  to  7  GPa  whereas 
Si02  demonstrates  this  behavior  at  >15  GPa  under  both  static  and  dynamic  compression 

[Kingma  et  al.,  1995].  As  shown  in  Figure  1,  transformation  of  a  dense  rutile-like  structure  (6.3 
g/cm3)  begins  at  ~4  GPa  and  is  complete  at  ~15  GPa.  The  shock  velocity-particle  velocity  plane, 
Figure  2,  explicitly  demonstrates  that  the  pressure  range  of  transformation  corresponds  to 
particle  velocities  in  the  300-800  m/s  range.  This  velocity  range  corresponds  to  particle 
velocities  imparted  by  small  arms  fire  projectiles  striking  GeC>2  glass. 

Earlier  motivated  by  the  substantial  increase  in  density  displayed  by  Hugoniot  curves, 
Chen  et  al.  [1999]  measured  the  release  isentropes  in  order  to  allow  calculation  of  the  excess 
energy  deposited  in  the  GeC>2  in  regions  where  the  vitreous  to  rutile  structure  phase  change 
occurred.  It  was  expected  that  this  hysteretic  phase  change  could  result  in  severe  attenuation  in 
the  amplitude  of  an  unsupported  diverging  shock  as  it  decays  through  the  14  to  4  GPa  range.  At 
lower  stresses  the  attenuation  is  expected  to  be  no  longer  anomalous.  In  silicate-igneous  and 
carbonate  and  other  sedimentary  rock  types,  the  attenuation  rate  when  fit  to  a  power  law  indicate 
the  shock  stress  cH  decays  with  radius  directly  beneath  a  point  source  at  the  surface  as: 

oH  l/r" 

where  n  varies  from  1.5  to  1.8  for  aH  ~  2-300  GPa.  [Ahrens  and  O'Keefe,  1985].  The  attenuation 
rate  from  fully  contained  explosions  over  a  similar  stress  range  is  similar  with  n  =  1.7  -  1.8 
[Cooper,  1977]. 

The  release  isentropes  (Figure  lb)  for  Ge02  glass  were  measured  using  the  experimental 
setup  of  Figure  4.  In  the  two  forward  experiments  (Figure  4a)  3-4  mm  thick  samples  were 
impacted  at  1.5  and  2.0  km/sec  with  2.6  mm  thick  tungsten  plates  (Table  1).  In  the  reverse 
experiments  (Figure  4b)  3  or  4  mm  thick  plates  of  Ge02  were  launched  at  0.8  to  1 .5  mm  thick 
6061  Al  buffer  plates  at  speeds  of  0.6  to  2.0  km/sec. 

A  typical  VISAR  record  is  shown  in  Figure  5a  (this  was  a  forward  shot).  The  wave 
profile  with  a  14  GPa  amplitude  displays  a  0.6  GPa  elastic  precursor  and  is  fit  with  a  hysteretic 
model  to  produce  Figure  5c.  Panel  5b  was  calculated  using  the  standard  Gruneiser  equation-of- 
state  but  lacking  a  hysteretic  equation  of  state.  The  peak  pressure  duration  that  is  calculated  is 
clearly  much  longer  than  observed  in  Figure  5a.  The  numerical  simulations  shown  were 
performed  with  the  one-dimensional  compressible  flow  computer  code,  WONDY  [Kipp  and 
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Lawrence,  1982].  We  compare  calculated  profiles  to  those  measured  with  the  VISAR  to  yield  a 
series  of  release  isentropes  that  are  shown  in  Figure  lb.  Only  the  initial  release  isentropes 
actually  are  well  defined  by  the  VISAR  data.  The  initial  unloading  is  induced  by  the  arrival  of  a 
rarefaction  wave  from  the  rear  of  the  flyer  plate  in  both  forward  and  reverse  configuration 
(Figure  4a, b).  The  particle  velocity  profiles  are  useable  until  arrival  of  lateral  rarefaction  waves 
at  the  VISAR  reflection  point  in  the  center  of  the  sample  about  0.5  fis  after  the  initial  rarefaction 
wave  from  the  rear  of  the  flyer  plate  attenuates  the  peak  shock  pressure. 

Using  the  release  isentropes  of  Figure  lb  for  the  phase  transition  region,  Chen  et  al.  [, 
1998  #3783]  calculates  the  stress  wave  attenuation  for  the  geometry  of  Figure  6.  We 
approximate  the  stress  attenuation  from  a  projectile  striking  a  target  by  examining  the  stress 
wave  decay  from  a  spherical  flyer  plate  impacting  a  sample  in  radial  geometry.  The  peak  stress 
wave  decrease  with  propagation  radius  is  shown  in  Figure  7  for  both  reversible  and  irreversible 
unloading  model.  It  is  clear  that  more  rapid  decay  of  peak  stress  amplitude  with  radius  takes 
place  in  the  pressure  regime  of  the  hysteretic  phase  change.  The  shock  pressure  from  the 
spherical  source  (which  is  a  model  of  the  stress  attenuation  along  the  centerline  beneath  a  normal 
projectile  impact  on  a  Ge02  half  space)  was  predicted  to  decay  in  Chen  et  al.  [,  1998  #3783]  with 
a  power-law  (Equation  1)  with  n  =  2.71 .  This  is  considerably  more  attenuative  than  the  regimes 
where  the  phase  change  does  not  occur.  We  point  out  that  analogous  behavior  over  the  shock 

stress  range  9-30  GPa  in  Si02  were  calculated  for  radial  flow  by  Swegle  [1990]  although  these 
results  were  not  fit  to  a  power-law  curve. 

PLANE-WAVE  PROPAGATION  IN  GeOz 

Using  Manganin  gauges,  in  planar  flow  geometry,  we  conducted  measurements  of  elastic 
precursor  decay  from  final  shock  states  in  the  range  of  8-16  GPa  using  the  geometry  of  Figure  9. 
Previously  Wackerle  [1962]  measured  elastic  precursor  decay  in  fused  quartz.  Because  the  fused 
quartz  elastic  Hugoniot  is  convex  upward  in  the  pressure- volume  phase  to  6  GPa  the  initial 
portion  of  the  Si02  glass  precursor  are  dispersive.  However,  at  stresses  above  the  6  GPa  level  a 
normal  steep  shock  front  is  induced  because  the  Hugoniot  becomes  concave  downward.  This 
occurs  up  to  the  Hugoniot  elastic  limit  (HEL)  of  9  GPa.  Thus  Si02  can  be  viewed  as  displaying 
a  partially  supported  elastic  precursor.  We  find  that,  as  expected,  this  wave  decay  occurs  slowly 
with  distance  and  goes  with  n  =  0.54.  In  contrast,  the  elastic  precursor  in  Ge02,  which  is  always 
dispersive  up  to  the  4  GPa  (Figure  10),  decays  more  rapidly  with  n  =  1.24  (Figure  1 1). 

SPHERICAL  WAVE  DECAY 

We  have  only  fired  two  experiments  that  unequivocally  demonstrated  unsupported  wave 
decay  in  spherical  geometry  using  the  set-up  shown  in  Figure  12.  Asa  result  of  the  radial  flow 
induced  by  the  spherical  impactor  the  manganin  gauge  resistance  (50  ohms-Dynasen-MN4-50) 
increases  as  a  result  of  both  stretching  as  a  result  of  the  radial  diverging  flow  and  a  result  of 
being  encompassed  by  the  increasing  shock  stress.  To  take  this  into  account,  we  employ  also  in 
the  same  planes,  Constantan  strain  gauges  (MN/CN-4-50-ER  Dynasen).  Constantan  alloy  is 


4 


relatively  insensitive  to  shock  heating.  However,  like  Manganin,  it  suffers  stretching  when 
placed  perpendicular  to  a  radial  flow. 

We  combined  data  from  Manganin  and  Constantan  and  employed  the  latter  to  apply  a 
stretching  correction  to  the  Manganin  record,  for  example,  Figure  13.  Spherical  elastic  precursor 
decay  for  one  experiment  is  shown  in  SiC>2  in  Figure  14  and  vitreous  GeC>2  in  Figure  15.  An 
approximate  fit  to  the  Ge(>2  data,  taking  into  account  a  region  0.2  X<>  in  radius  (X0  is  the 
spherical  impactor  radius  of  Figure  11).  In  this  close-in  (region  A),  Figure  15, 0<x/Xo<0.2  the 
stress  wave  attenuation  is  slight  as  the  flow  is  nearly  planar  and  only  a  slight  deviance  of  stress 
occur  with  increasing  radius  (due  to  the  slightly  diverging  flow).  This  feature  of  shock 
attenuation  from  impact  was  observed  by  Ahrens  and  O’Keefe  [1977]  (Figure  16). 

In  Region  B  (Figure  15)  the  shock  attenuation  is  dominated  by  the  hysteretic  phase 
transition  and  a  decay  rate  of  N=3.35  is  estimated  in  Figure  17. 

CONCLUSIONS 

The  stress  wave  attenuation  from  spherical  unsupported  shocks  in  engineering  and 
geologic  materials  yield  stress  wave  decay  exponents  of  n  =  1.5  to  1.8  (Equation  1)  in  the  0.5  to 
150  GPa  pressure  range.  The  present  data  are  consistent  with  hysteretic  behavior  associated  with 
n~3.  We  calculate  a  value  of  n  =  3.35  on  the  basis  of  two  spherical  experiments  within  the  phase 
transition  region  with  virtually  no  control  of  uncertainty.  This  compares  to  the  n  =  2.7  predicted 
exponent  [Chen,  1998  #3783]  in  this  stress  range.  GeC>2  appears  to  be  representative  of  a 
hysteretic  phase  change  materials  and  would  appear  to  be  extremely  attenuative  in  the  4-15  GPa 
phase  change  regime.  Further  work  on  a  longer  time  scales  and  their  effect  on  ballistic 
penetration  must  still  be  studied. 
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TABLES  AND  FIGURES 


Table  1 

VISAR  experiments  on  GeC>2 
1.  Forward  shots 


Flyer  plate  (Tungsten) 

Stationary  plate  (GeCh) 

Buffer  (A 16061) 

Shot 

Thickness 

Velocity 

Thickness 

Density 

Thickness 

No. 

(mm) 

(km/s) 

(mm) 

(Mg/m3) 

(mm) 

953 

2.627 

1.53 

3.835 

3.632 

0.743 

±0.003 

±0.04 

±0.002 

±0.002 

±0.006 

955 

2.645 

1.99 

4.059 

3.631 

0.923 

±0.001 

±0.05 

±0.002 

±0.003 

±0.001 

2.  Reverse  shots 


Flyer  plate  (GeC^) 

Driver  plate  (A 16061) 

Shot 

Thickness 

Velocity 

Density 

Thickness 

No. 

(mm) 

(km/s) 

(Mg/m3) 

(mm) 

957 

4.087 

1.47 

3.630 

1.459 

±0.001 

±0.05 

±0.04 

±0.001 

958 

4.128 

1.96 

3.633 

1.354 

±0.001 

±0.08 

±0.03 

±0.001 

965 

3.029 

0.62 

3.618 

0.774 

±0.001 

±0.01 

±0.03 

±0.001 

Densities  of  tungsten  and  A 16061  are  19.2  and  2.70  Mg/m  ,  respectively 
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Figure  1.  Hugoniots  and  release  isentropes  for:  (a)  reversible  phase  change  models;  (b)  irreversible 
(several  measured  releases  at  different  pressures  are  shown).  For  reversible  model,  release  are  not 
distinguishable  from  Hugoniot.  “Waste  heat”  is  area  between  the  Rayleigh  line  (dashed)  and  the  solid 
curve  release  isentrope  induces  normal  attenuation.  The  hysteretic  behavior  of  Hugoniot-release 
isentrope  loop  in  (b)  enhances  dissipation  of  shock  energy  and  gives  rise  to  increased  attenuation. 


Figure  2.  Shock  pressure  versus  volume  of  vitreous  Ge02.  SL  and  HL  are  Hugoniot  states  shock 
compression  and  isothermal  (hydrostatic)  compression  states,  respectively. 
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Particle  velocity  (km/s) 

Figure  3.  Shock  wave  versus  particle  velocity,  vitreous  Ge02,  initial  density. 


(a)  Forward  Experiment 


(b)  Reverse  Experiment 


Figure  4.  Diagram  of  Ge02  release  isentrope  experiments:  (a)  Metal  flyer  plate  impacts  sample  assembly 
and  generates  a  shock  wave  through  the  sample.  The  rearward  propagating  shock  wave  upon  passing 
through  the  flyer  plate  reflects  the  flyer  foam  interface  and  propagates  forward  as  a  rarefaction  wave. 
VISAR  probe  detects  particle  velocity  change  at  aluminum  reflector-window  interface  due  to  shock  and 
following  rarefaction  waves,  (b)  Reverse  experiment  is  similar,  except  the  flyer  plate  is  the  sample  under 
study  (Ge02),  and  it  directly  impacts  the  buffer/reflector.  Reverse  experiments  usually  provide 
continuous  release  data  to  lower  pressures,  but  Ge02  is  prone  to  damage  during  projectile  launch. 
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Figure  5.  Particle  velocity  versus  time  measured  by  VISAR  (a)  and  calculated  (b,c);  (b)  WONDY  fit 
experiment  using  a  reversible  phase  change  model;  (c)  WONDY  fit  to  experiment  using  an  irreversible 
phase  change  model. 
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Figure  6.  Configuration  for  calculating  shock  wave  decay  in  spherical  geometry.  At  time  zero,  tungsten 
core  impacts  surrounding  Ge02  at  1 .5  hr/sec  simulating  impact  at  surface  of  a  half-space. 


Figure  7.  Shock-wave  profiles  at  a  series  of  times  (separated  by  1  p  sec)  in  spherical  symmetry:  (a) 
Reversible  model;  (b)  Irreversible  model.  In  geometiy  of  Figure  6.  Impactor:  tungsten,  radius  from  6  to 
8  mm,  moving  at  1.5  km/s;  target:  Ge02  (reversible  model),  24  mm  thick. 
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Figure  8.  Logio  shock-wave  amplitude  versus  logio  radius  (spherical  symmetry).  Results  for  reversible 
and  irreversible  models  are  shown,  each  with  impact  velocities  of  0.5  km/s  and  1.5  km/s.  Only  1.5  km/s, 
with  hysteretic  behavior  results  in  excess  stress  wave  attenuation  over  15  to  4  GPa  range. 


Figure  9.  Configuration  for  propagation  experiments  in  one-dimensional  planar  flow  in  the  Ge02  glass. 
Flyer  plate  material  is  2024  A1  or  Cu  with  thickness  4  to  6  mm.  Manganin  stress  gauges  are  from 
Dynasen,  Inc.  Ramp  elastic,  deformation  and  also  release  states  are  calculated  from  shock  particle  and 
shock  propagation  velocity  and  Riemann  integral  relation.  a(t)  =  A0  +  AidR(t)  +  A2dR2(t)  +A3dR3(t)  + 
AtdR4(t)  [Rosenberg  et  al.,  1980].  Here  dR(t)=R(t)-Ro)/Ro(R(t)  is  gauge  resistivity  during  loading  at  time 
t  and  Ro  is  initial  gauge  resistivity.  Pressure  is  calculated  from  measured  voltage  change  using  calibration 

constant  Aj.  When  c  <  1.5  GPa,  Ai=50  (GPa)  and  other  coefficients  are  zero.  When  1.5  GPa,  A0  = 
0.572  GPa,  A,=29.59  GPa/ohm-m,  A2=95.20  GPa/(ohm-m)2,  A3=3 12.74  GPa/(ohm-m)3  and  A4-331.77 
GPa/(ohm-m)4. 
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Figure  10.  Plane-shock  and  release  wave  profiles  in  vitreous  GeC>2  development  of  4  GPa  elastic 
precursor  with  ramp  structure  is  consistant  with  convex  upward  elastic  Hugoniot. 
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Figure  11.  Hugoniot  elastic  shock  amplitude  versus  normalized  propagation  distance  for  vitreous  SiC>2 
and  GeC>2.  Normalization  distance  is  initial  thickness  of  flyer  plate. 


Figure  12.  Diagram  of  spherical  impactor  radial  wave  attenuation  experiment.  Each  gauge  position 
incorporates  both  a  Manganin  gauge  (sensitive  to  shock  pressure  and  stretching  due  to  radial  flow)  and  a 
Constantan  strain  gauge  which  measures  stretching  arising  from  radial  flow. 
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Figure  13.  Stress  wave  profiles  in  spherically  diverging  flow.  Manganin  stress  gauge  detects  increasing 
resistance  with  time.  Signal  is  superposition  of  increased  resistance  from  increase  in  pressure  of  shock 
plus  stretch  of  manganin  gauge  produced  by  (lateral)  strain  parallel  to  spherical  wave  front.  Lateral  strain 
is  measured  using  a  Constantin  strain  gauge.  The  Constantin  strain  gauge  record  allows  correction  of 
stress  gauge  (Manganin)  signal  so  as  to  take  into  account  stretching. 


Figure  14.  Stress-wave  profiles  for  spherical  impact  in  fused  SiC>2.  Peak  stress  measured  is 
approximately  at  amplitude  of  Hugoniot  elastic  limit.  Hence  observed  data  give  a  measure  of 
unsupported  elastic  shock  attenuation  in  fused  SiC>2. 
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Time  (jxs) 


Figure  15.  Stress-wave  profiles  for  spherical  impactor  in  vitreous  GeC>2.  Arrows  indicate  the  amplitude 
of  the  precursor  originating  from  impact  surface.  Note  the  peak  stress  decays  for  8  to  4  GPa  over  a  3  mm 
increment  in  radius. 


Figure  16.  Logio  (peak  shock  presure)  versus  logio  (normalized  radius)  at  various  impact  velocities  for 
gabbroic  anorthosite  impactors.  Arrows  indicate  one-dimensional  flow  pressures.  Close  in  region  that 
has  minor  attenuation  results  from  a  quasi-planar  flow  (Regime  A  of  Figure  17).  Minor  divergence  of 
flow  induces  minor  stress  wave  decay  immediately  beneath  striking  projectile  (after  Ahrens  &  O’Keefe 
1977) 
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Figure  17.  Peak  shock  pressure  versus  propagation  distance  A:  near  field  attenuation.  B:  phase 
transition  dominated  decay.  Here  peak  stress  decay  is  1/r"  where  n  =  3.35.  C:  ramp-wave  dominated 
decay.  D:  elastic-wave  decay. 
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Abstract 

This  thesis  applies  recently  developed  continuum  theories  of  diffusionless  phase 
transformations  in  solids  to  the  study  of  impact  problems  involving  materials  which 
can  experience  such  phase  changes.  Our  objective  is  to  compare  the  theoretical 
predictions  against  certain  experimental  results. 

In  the  experiments  of  interest,  a  face-to-face  impact  occurs  between  a  disk  of 
amorphous  germanium  dioxide  and  another  material,  either  tungsten  or  an  alu¬ 
minum  alloy.  The  Ge02  is  believed  to  transform  to  another  phase  if  sufficient 
compressive  stress  is  achieved. 

We  model  these  experiments  using  one-dimensional  finite  elasticity.  Phase¬ 
changing  materials  are  represented  by  non-convex  potential  energy  functions.  This 
can  produce  phase  boundaries  that  propagate  subsonically  or  supersonically  with 
respect  to  the  slower  longitudinal  wave  speed  of  the  two  phases.  When  a  subsonic 
phase  boundary  is  possible,  it  is  not  uniquely  determined  by  the  fundamental  field 
equations  and  jump  conditions.  Uniqueness  is  obtained  by  invoking  a  nucleation 
criterion  to  control  the  initiation  of  the  new  phase,  and  a  kinetic  relation  to  govern 


V 


its  evolution. 

The  experiments  considered  here  axe  sufficiently  long  in  duration  («  3  ps)  that 
several  reflections  and  wave  interactions  occur,  and  the  analysis  becomes  analytically 
intractable.  Accordingly,  a  finite-difference  method  of  Godunov  type  is  employed 
to  analyze  these  experiments  numerically.  Methods  of  Godunov  type  treat  adjoin¬ 
ing  discretized  spatial  elements  as  the  two  sides  of  a  Riemann  problem,  which  is 
typically  solved  approximately  by  linearizing  around  the  initial  conditions  on  each 
side.  Fortuitously,  all  constitutive  models  employed  in  this  thesis  are  such  that  the 
required  Riemann  problems  can  be  solved  exactly  without  too  much  effort. 

Simulations  utilizing  the  numerical  method  demonstrate  that  the  impact  re¬ 
sponse  of  a  material  is  sensitive  to  the  kinetic  relation  that  enters  the  model.  It 
appears  the  theory  may  offer  a  plausible  description  of  the  experiments,  though  the 
restrictions  placed  on  the  constitutive  models  herein  seem  too  severe  to  provide  a 
good  quantitative  match  to  the  experimental  results. 
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Abstract 

Stress  wave  profiles  in  vitreous  GeC>2  were  measured  using  piezoresistance  gauges  in 
the  pressure  range  of  5  to  18  GPa  under  planar  plate  and  spherical  projectile  impact. 
Experimental  data  show  that  the  response  of  vitreous  GeC>2  to  planar  shock  loading 
can  be  divided  into  three  stages:  (1)  A  ramp  elastic  precursor  has  peak  amplitude 
of  4  GPa  and  peak  particle  velocity  of  333  m/s.  Wave  velocity  decreases  from  initial 
longitudinal  elastic  wave  velocity  of  3.5  km/s  to  2.9  km/s  at  4  GPa;  (2)  A  ramp  wave 
with  amplitude  of  2.11  GPa  follows  the  precursor  when  peak  loading  pressure  is  8.4 
GPa.  Wave  velocity  drops  to  the  value  below  bulk  wave  velocity  in  this  stage;  (3)  A 
shock  wave  achieving  final  shock  state  forms  when  peak  pressure  is  >  6  GPa.  The 
Hugoniot  relation  is  D  =  0.917  +  1.711u  (km/s)  using  present  data  and  the  data  of 
(  Jackson  and  Ahrens  [1979]  when  shock  wave  pressure  is  between  6  and  40  GPa  for 

Po  =  3.655  g/cm3.  Based  on  the  present  data,  the  phase  change  from  4-fold  to  6-fold 
coordination  of  Ge+4  with  O-2  in  vitreous  GeC>2  occurs  in  the  pressure  range  of  4  to 
15  ±  1  GPa  under  planar  shock  loading.  Comparison  of  the  shock  loading  data  for 
fused  SiC>2  to  that  on  vitreous  GeC>2  demonstrates  that  transformation  to  the  rutile 
structure  in  both  media  are  similar.  The  Hugoniots  of  vitreous  GeC>2  and  fused  Si02 
are  found  to  coincide  approximately  if  pressure  in  fused  SiC>2  is  scaled  by  the  ratio 
of  fused  Si02  to  vitreous  GeC>2  density.  This  result,  as  well  as  the  same  structure, 
provides  the  basis  for  considering  vitreous  GeC>2  as  an  analogous  material  to  fused 
SiC>2  under  shock  loading.  Experimental  results  from  the  spherical  projectile  impact 
demonstrate:  (1)  The  supported  elastic  shock  in  fused  SiC>2  decays  less  rapidly  than 
a  linear  elastic  wave  when  elastic  wave  stress  amplitude  is  higher  than  4  GPa.  The 
supported  elastic  shock  in  vitreous  GeC>2  decays  faster  than  a  linear  elastic  wave;  (2) 
In  vitreous  GeC>2,  unsupported  shock  waves  decays  with  peak  pressure  in  the  phase 
/  transition  range  (4  -15  GPa)  with  propagation  distance,  x,  as  oc  1/x-3-35,  close  to 
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the  prediction  of  Chen  et  al.  [1998].  Based  on  a  simple  analysis  on  spherical  wave 
propagation,  we  find  that  the  different  decay  rates  of  a  spherical  elastic  wave  in  fused 
SiC>2  and  vitreous  GeC>2  is  predictable  on  the  base  of  the  compressibility  variation 
with  stress  under  one-dimensional  strain  condition  in  the  two  materials. 


